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Abstract 

The surface stress can have important effects on the elastic properties 
of nano-sized structures. Here we analyze the effect of surface stress 
on Mode-3 crack displacement and stress-field solutions under the as- 
sumption of linear elasticity. We show that surface effects generate 
non-K terms near the crack-tip. We also find the effect of such terms 
on the conventional small-scale yielding assumption. 

1 Introduction 

The distinction between the energy involved in creating a new surface, as 
compared to stretching an already existing one, has been reemphasized many 
times over the years 121 El H] • While the former can fairly be taken as 
a constant per unit area, the latter may have complicated dependence on 
surface strain. Here we try to find the effect of the latter ('surface elasticity') 
on crack field solutions when this dependence is quadratic akin to its bulk 
counterpart. 

Physically, this effect will be important when the solid dimensions are 
of the order of k/h where k and \x are surface and bulk elasticity moduli re- 
spectively. For real materials this ratio is of the order of nanometers. Since 
the dimensions of cracks in many materials such as glass range from few 
angstroms to micrometers, and crack elastic field solutions depend largely 
on surface tractions, it seems worthwhile to analyze the effect of surface 
elasticity on crack field solutions. At the outset, the use of continuum elas- 
ticity at such dimensions may be questioned. Some years ago Shenoy et 




Figure 1: The specimen 



al compared the continuum elasticity results with atomistic simulations |Hj 
and were able to show that bending and torsional behavior of a solid at such 
dimensions can be accounted for by taking surface elasticity into account. 
Here in similar spirit we try to formulate the problem of a mode-3 crack for 
a material having elastic surfaces. 

2 Elastic solution for a Mode-3 Crack with Surface 
Elasticity 

In this section we solve for the elastic field of a mode-3 crack for a solid 
having an elastic surface. We follow a perturbative approach with the field 
in the absence of surface elasticity as the unperturbed solution. The domain 
0, consists of < r < R, —ir < 9 < tt (see Fig 1). Let the loading at the 
external boundary r = R be of the form T rz (R,9) = T a F{9). 

In the absence of surface tension and surface elasticity effects, the only 
non-zero displacement w (r,9) in the z-direction satisfies 

V 2 w (r,9) = (0<r<R) (1) 
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~ = (r,e = n,-n) (2) 
r ot) 

^ = r a F(6) (r = R) (3) 

Supplementing the above equations with the condition of finite displace- 
ment at r = 0, the anti-symmetric solution (i.e. the one which contributes 
to crack-opening) to above set of eqns. can be easily found as 

oo 

w o (r,0) = ^2 w on (4) 

n=0 

w on = A n r n+l / 2 sin{(n + 1/2)6) (5) 

A - TaFn (ft) 

H(n + 1/2)^-1/2 W 

where F n = I F(6)sin((n + 1/2)6) d6. 

A general formulation of elastic surfaces was given by Gurtin and Mur- 
doch jSj. Following (Sj, in the presence of surface elasticity /tension effects, 
the equation of equilibrium in the bulk remain unchanged (i.e. same as 
Eq. l[T|)) while the boundary conditions at a free surface are changed from 
(JijUj = to: 



S a /3 = Po<>a,P + (A + Cr )u-y^5 a: f3 + (fl - (J ){u a ^ + Uj3 )0l ) + &oU a ,f3 

Here S a p denote surface stress tensor while A, fj, a and o Q are surface 
elastic moduli. If the normal to the surface lies along direction 2, then Latin 
subscripts i, j, k range over 1,2,3 while Greek subscripts a, (5 range over only 
the directions lying along the surface viz. 1 and 3. The derivatives such 
as Sia t a or u a ^ employed in the above eqns. are surface derivatives i.e. 
projection of usual derivatives along the surface. Note that S is linear in 
surface strain E s where E s a p = \{u a ^ + up >a ). 

We switch to cartesian coordinates for a while for sake of mathematical 
convenience. Along the crack faces (6 = tt, — ir), the negative x axis coincides 
with 6 = 7r(— tt) for y = 0+((T). 

Substituting u a = w(x,y)5 a 3 in above equations 

1 dw , 

E s = - ir ei®e 3 + e 3 §ei (7) 

2 ox 

dw dw 
S = cr G ei 0ei + - ct ) — (ei (g>e 3 + e 3 (giei) + CT — ei (g)e 3 (8) 
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Taking surface divergence of S one obtains 

d 2 w 

■) (9) 

where k = fi — a . The first term arises because the constant part of surface 
stress in Eq. © drops abruptly to zero at the crack tip. Its effect on stress 
fields was solved for many years ago by Thomson et al |0j who showed that 
such terms lead to 1/r singularity in stresses at the crack-tip. As we are 
concerned with the effect of surface elasticity and not surface tension in this 
work and since the overall field can be found as superposition of the fields 
due to both terms individually, we drop the first term in above eqn. Keeping 
this in mind, the above eqn. in r — 9 coordinates may be written as 

1 dw k d 2 w q 
r 89 n dr 2 
Idw nd 2 w 

-— - = 7T~n~ " = (11) 

r o9 \i or 1 

Eqn. Q in conjunction with (|10j) and 1)11(1 completes the specification 
of boundary value problem for our problem. We were unable to find a 
straightforward method to solve this set of equations. We will solve for w as 
a series expansion in terms of the parameter . We make the assumption 

A « L 

Next write w(r, 9) as 



w{r,9) = ^w n (12) 

ra=0 

w n = w on + (K/fi)gin + (K/n) 2 g 2n + ■■■ + (K/n) k g kn + ■■■ (13) 



where w on is given by (JSJ). Throughout our analysis we assume that 
Since u> n 's are independent, each must satisfy the equation V 2 w n = in the 
bulk and the boundary conditions l)l()j). (fTT)) individually. Thus one gets 
the following set of equations 



V 2 2fcnM) = (0<r<R) 

for all k, and 

1 dw on 

— = 

r 89 
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1 dgin 
r 39 
1 dg 2n 
r 89 



d 2 w, 



on 



(14) 



dr 2 




dr 2 



71 



{0 = vr,-7r} 



(15) 



1 dgu 
r 89 



O 2 



Qk-Vn 
dr 2 



where upper sign stands for 9 = tt and lower for 9 = —tt. Substituting 
for w on from © in (|14|) and solving for g\ n 



9ln = ( n + l/2)A n r n - 1 l 2 cos{{n-l/2)9) + a n r n - 1 l 2 sin{{n-l/2)9) 



where a n is an arbitrary coefficient. Since the loading due to surface elas- 
ticity is symmetric (note the sign change in boundary condition at 9 = tt 
and —tt) , all the terms containing k should be symmetric in 9. Therefore 
a n = 0. Thus the terms containing k do not contribute to the crack opening 
displacement but the stress field is necessarily altered due to surface elastic- 



Again by the same argument g 2n is seen to be identically equal to zero. It 
immediately follows that g 2n = <?3n = ... = 0. Thus the perturbative solution 
is 



Wn = A n r n+l ' 2 sin{{n + 1/2)9) + {k/h){ti + l/2)A n r n -^ 2 cos({n - 1/2)9) (16) 

3 Small Scale Yielding and Surface Elasticity 



From equations (|12|) , (|16|) one observes that the displacement diverges if we 
include the term n = 0. As and Ruina [7] have shown through their anal- 
ysis that one needs to be careful while discarding non-iT singular(and non- 
singular) terms for crack problems. In fact, such terms are always present 
if one take a more realistic situation where there is a nonlinear zone £1' sur- 
rounding the crack as shown in figure 2. Here we do an analysis similar 
to [7j to determine the condition for small-scale yielding in the presence of 
surface effects. 

Assume that the nonelastic zone is centered at the origin x, y = and has 
dimensions smaller than the crack length R. We now consider a completely 



ity 

Substituting the expression for g\ n in 



G5J 



1 dg 2 , 
r 80 



•n 



= 
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Figure 2: The specimen 



elastic problem in the domain p < r < R, p being the outermost radius of 
0,'. The traction at the inner radius r = p is taken as 

Trz(p,0) = r H e 

together with the loading r a F(6) at r = R. The antisymmetric solution to 
the elastic problem in the absence of surface effects can now be found as 

oo 

w (r,6) = ^2 w on 

n=0 

Won = {C n r n+1 l 2 + D n r- ( - n+1 ^}sin{(n + 1/2)6) 

where F n = I f^F(0)sin((n + 1/2)9) d6 and e = p/R. 

Following the notation of 0, the above field can be expressed as a sum 
of two terms: w AI and w AO where AI stands for Antisymmetric, Inner 
traction free and AO for Antisymmetric, Outer traction free. For sake of 
completeness we write the expressions for these fields: 

oo 

Wo (r,e) = + 

n=0 
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ul r n+l/2 n 2n+1 

vAi = 11 "J - -7^ + ^-—^W(n + 1/2)0) 

~ ifl ~.n+l/2 pn+3/2 
AO T °n r n+3/2 r ' , n+3/2 n 



/x(n + l/2) L fl n -V2 r n+l/2J vv / ^ 



with e = p/R, 
and 

6* = F n /[1 - e 2n+1 ] with F n = i F(0)ain((n + 1/2)9) d6, 
and 

ft 2 = -H n /[1 - e 2n+1 ] with H n = \ jZ n H(9)sin((n + 1/2)6) d6, 

In the presence of surface effects the boundary conditions (J2J), (J2J) remain 

unchanged except that now our domain is p < r < R, 6 = tt, —it. Following 

exactly the same strategy as before, we obtain the following perturbative 

solution: 



w(r,6) = ^2,w n 

n=0 



w„ 



T r n-l/2 n 2n+1 

< + <° + ^(K/^)b 1 n [-— Tr2 cos((n - 1/2)6) + J ] x 



V"/H"n[jj,.i/2"™U" i ftn-l/2 r n+3/2> 

2 r.n+3/2. 



T r n-l/2 

cos((n + 3/2)0) + -Z(k/M [ e ™+3/2___ cos((ri _ 1/2)0) + 

pn+3/2 

e^-^]cos((n + 3/2)6) 

From this, the stress fields are readily obtained, 
dw 

Trz = 

or 

oo n— 1/2 2n+l 

= EfK^I^ - ^-V2 r n + 3/2 ]^((^ + 1 /2)^)} + {r («//i)6i x 

n=0 

r n-3/2 n 2n+1 

[(» " 1 /2) ] ^c OS ((n - 1/2)0) - (n + 3/2) _^_^ ]co8 (( n + 3/2)0)} 

n-1/2 , on+3/2 

+ {r # [ e ™+3/2___ _ e n + 3/2__ ]s in((n + 1/2)9)} + {To{K/fl)b 2 n x 

r n-3/2 on+3/2 

[(n - 1/2)^^-^008^ - 1/2)6) - (n + 3 /2) e " +3 / 2 -^] COS ((n + 3/2)0)}] 
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and 



1 dw 

oo r n-l/2 p 2n+1 

= E[Wn[jp=I72 + fi n-l/2 r n + 3/2 ] COS ((" + W)} ~ X 
n=0 

r n~3/2 „2n+l 

[(» " 1 /2)^T 7 ^^n((n - 1/2)0) + (n + 3/2)^^]™((n + 3/2)0)} 

-n-1/2 pn+3/2 
+ {r ^[ e n+3/2___ + e n + 3/2£L__ ]cos((n + 1/2) ^ )} _ {^/^ x 

r n-3/2 dti+3/2 

[(n - l/2)e n+3 / 2 j^sin((n - 1/2)9) + (n + 3/2)e™+ 3 / 2 -^^] S m((n + 3/2)0)}] 



From above equation the condition for K term) field dominance can 
be determined. Comparing the total stress field due to K field to that of 
next singular term one obtains the following expression for the region 
in which K field dominates: 



[(T a /T )(p + a ) + e 3 / 2 (R + a )] 
[( Ta /r )(l + a /R) + eV 2 (l+a e/R)) 1 ) 

where a = n/fi 

It immediately follows from above eqn. that when r Q = 0, the K-dominance 
region diminishes in its radius by a Q while for the case when r a = the 
-K= term always dominates the K term as it does in the absence of surface 

ryr J 

elasticity. 



4 Discussion 



To put our analysis in perspective, many articles pertaining to effect of 
surface stress on crack field solutions have appeared in the past H3 . 
The main difference with our work is that we have considered surface stress S 
as having linear (and thus surface energy quadratic) dependence on surface 
strain E s instead of being independent of it. As mentioned earlier that a 
constant surface stress leads to r" 1 singularity for a Mode-3 crack [Sj. We 
have shown that the linear term causes a r~ 3 / 2 singularity at the crack-tip 
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for e = 0. Physically speaking, non-i^T singular terms due to surface elasticity 
should be dominant only in the region of extent a Q = n/fj,. Indeed we see 
from Eq. (|17|) that surface elasticity decreases the extent of ET-dominant 
region by a Q . 

Though we are not aware of any direct experimental evidence which fa- 
vors a quadratic dependence of U s on E s over a linear one, Miller et al [3] 
showed that a quadratic dependence was consistent with atomistic simu- 
lations for simple problems of elongation, bending and torsion. Note that 
the variable S, which is responsible for size-dependent elastic properties in 

has same significance as k in this article under Gurtin-Murdoch's for- 
malism Finally, the Green's function obtained by Thomson et al [2] for 
symmetric dipole acting along a Mode-3 crack's surface may also be used 
for obtaining the results in this article but since in our problem the loading 
itself is dependent on the displacement field, we found it more expedient 
to directly solve the equations of elasticity instead of a solving a singular 
integro-differential equation. 
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